Tricritical behavior in itinerant quantum ferromagnets 
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It is shown that the peculiar features observed in the low-temperature phase diagrams of ZrZn2, 
UGe2, and MnSi can be understood in terms of a simple mean-field theory. The nature of the 
ferromagnetic transition changes from second order to first order at a tricritical point, and in a 
small external magnetic field surfaces of first-order transitions emerge which terminate in quantum 
critical points. This field dependence of the phase diagram follows directly from the existence of the 
tricritical point. The quantum critical behavior in a nonzero field is calculated exactly. 

PACS numbers: 



The ferromagnetic transition at the Curie points of 
the elements Fe, Ni, and Co is one of the best-known 
examples of a second-order phase transition. It is well 
understood in terms of the band theory of metals in con- 
junction with the theory of phase transitions Q. Recent 
experimental studies of ferromagnetic compounds with 
much lower Curie temperatures, among them ZrZn2 0, 
UGe 2 |3, andMnSi HQ, show enigmatic behavior which 
does not seem to fit into this well-established picture: If 
the low Curie temperature is further decreased by means 
of pressure-tuning, the nature of the transition changes 
from second order to first order at a tricritical point, and 
in a small external magnetic field surfaces or "wings" 
of first-order transitions emerge which extend from the 
coexistence line at zero field and terminate in quantum 
critical points. These regions of first-order transitions 
end in lines of critical points which are reminiscent of 
conventional liquid-gas critical points, and which connect 
the tricritical point with a quantum critical points in the 
zero-temperature plane. These observations are summa- 
rized in the schematic phase diagram shown in Fig. ^ 

This structure of the phase diagram is very remark- 
able, for two reasons. First, the ferromagnetic transition 
in zero field at high-temperature Curie points, most no- 
tably in the elemental ferromagnets, is invariably of sec- 
ond order. Also, Hertz's theory of the quantum ferromag- 
netic transition at T = [(| , and its extension to nonzero 
temperature Q, predict the ferromagnetic transition to 
be generically of second order. Second, the persistence 
of the first-order transition away from the zero-field axis, 
and the existence of quantum critical points at h ^ 0, 
came as a surprise . Yet the observed structure of the 
phase diagram as sketched in Fig. ^ seems to be generic, 
as demonstrated by the case of ZrZn 2 as the latest ex- 
ample, where a tricritical point emerged once sufficiently 
clean samples were produced 0- 

In this Letter we show that all of these observations 
can be explained by a theory that takes into account the 
fact that, in metallic systems at low temperatures, the 
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FIG. 1: Schematic phase diagram in the temperaturc- 
pressure-magnetic field (T—p — h) space. Shown are the ferro- 
magnetic (FM, dark shaded) and paramagnetic (PM) phases 
at h = 0, the tricritical point (TCP), and the two quan- 
tum critical points (QCP). Also shown are various lines of 
first (dashed lines) and second-order (solid lines) phase tran- 
sitions, and the "wing" surface of first order transitions (light 
shaded) . 



particle-hole excitations characteristic of systems with a 
Fermi surface couple to the fluctuations of the magnetic 
order parameter and substantially change the nature of 
the phase transition compared to the conventional theory 
ISLIlTJ. Furthermore, we identify the universality classes 
for all finite-temperature critical points in the phase di- 
agram, and we determine the exact critical behavior at 
the quantum critical points. 

Within the framework of this theory, a mean-field the- 
ory for three-dimensional systems is defined by a free 
energy density 

f=-fuf> + t<j) 2 + vcj) 4 In (0 2 /m 2 + T 2 /T 2 ) + uaS 4 (1) 

in terms of an order parameter 4>. Here t and it, v > 
are coefficients analogous to those in Landau theory • 
For later reference we note that within Stoner theory 
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[T3. , which is a particular realization of Landau the- 
ory, t = 1 — r t iVp, with Afa the density of states at the 
Fermi surface and T t a microscopic spin-triplet interac- 
tion amplitude, u is proportional to the second derivative 
of the density of states, and v — 0. too is a microscopic 
magnetization (e.g., one Bohr magneton /^b per volume 
of a unit cell), and To is a microscopic temperature (e.g., 
the Fermi temperature). The physical value of 4>, which 
we will denote by ip, is the one that minimizes /. For 
h = 0, tp is equal to the magnetization to, and for this 
case the free energy given by Eq. JJ| was first considered 
in Ref. ITil For h > 0, a derivation along the same lines 
as in Refs. 0, ^3 shows that the magnetization is now 
related to the physical value of the order parameter by 



m = ip- (4/4/r t )h. 



(2) 



The equation of state, which relates to, T, and h, is ob- 
tained by minimizing /, which leads to 



h = 2t(p + Av(p 3 In {if 2 / ml + T 2 /T 2 ) 

ip 2 /m 2 



+v v? 



ip 2 /m 2 +T 2 /T 2 



Amp 3 (3) 



in conjunction with Eq. (J2J. 

Notice that, at T = 0, both the free energy and the 
equation of state are nonanalytic functions of the order 
parameter by virtue of the logarithmic term. This is 
in sharp contrast to ordinary Landau theory , where 
/ is an analytic function of 0, and it reflects the fact 
that the particle-hole excitations have been integrated 
out to obtain a free energy in terms of the order param- 
eter only. This is a particular example of a more general 
phenomenon, see Ref. 0. 

In order to discuss the mean-field theo ry g iven by Eq. 
(PJ, we first recall the solution at h = There is a 
tricritical point at (t = 0, T = T tc ), with T tc = T e~ u / 2v . 
At T — 0, the transition occurs at t = t\ = m 2 e~ 1 v e~ u l v 
and is of first order with the magnetization changing dis- 
continuously from to = to to = m\ = too e~ 1 l 2 e~ u l 2v . 
The line of first-order transitions at temperatures T < 
T tc is determined by f(tp) — f'{tp) — and can be given 
explicitly in parametric form (0 < s < 1): 



t(s) = ii(l-s) 2 
T(s) 



T tcS 1 /2 e -d- S )/2_ 



(4a) 
(4b) 



We now turn to the properties of the equation of state 
for h ^ 0. Consider first the T — plane. Suppose 
t has been tuned to t\, so that at h = the free en- 
ergy has two degenerate minima, one at (p = and one 
at ip = mi. For a small h > this double-minimum 
structure persists, and the two minima can still be made 
degenerate by increasing t, see Fig. There thus still 
is a first-order transition. However, with increasing h 
the two minima merge at a point where the first three 
derivatives of / vanish: f'((p c ) — fife) — f'ifc) = 




FIG. 2: The free energy at T = for h = 0, t = ti (solid 
curve) and h = 0.02, t = 0.1878 (dashed curve), respectively. 
In both cases, u — v — mo = 1 and T — 0. 
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FIG. 3: The surface of first-order transitions in the space 
spanned by T, t, and h. It is bounded by lines of first-order 
transitions in the T — and h — planes, respectively, and 
by the line of second-order transitions discussed after Eqs. 
A symmetric surface extends into the region where h < 0. 



with ip c 



toq e 



-13/12 e u/2v_ This condition determines 



a critical point (t c , h c ) in the T = plane that marks the 
endpoint of a line of first-order transitions. One finds 



t c = 6v(pl = 6e- 13/6 > 



,-u/v 



16 3 
= ~J Vifc 



16 



-13/4, 



-3u/2v 



(5a) 
(5b) 



This discussion can be repeated for any fixed value of 
T < T tc . Accordingly, there is a line of critical points 
connecting the tricritical point at (T = T tc , t = 0, h = 0) 
and the quantum critical point at (T = 0, t = t c , h = h c ). 
A parametric representation for this line is 

t(s) = ti(4s 2 + 5s + 6)(l-s) 2 e 1+ff(s) , (6a) 
T(s) = T tc s 1/2 e ff(s)/2 , (6b) 



h(s) = +* + !)(! -*) ' e 



5/2 13/4+3 ff (s)/2 



(6c) 



where 



9(s) = — (4 S 2 + 7 S + 13)(1- S ). 



(6d) 



This line forms a boundary of a surface of first-order tran- 
sitions that is shown in Fig. |3| 
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We see that the phase diagram obtained from Eq. Q 
has the same structure as the one observed experimen- 
tally, see Fig. As is the case with Landau theory, the 
phase diagram in a space spanned by observables will be 
a stretched and rotated version of the one in T — t — h 
space, since the parameters of the theory are complicated 
functions of pressure, temperature, and magnetic field 
However, a quantitative comparison with experi- 
ment can be made by expressing the zero-temperature 
critical field strength h C: Eq. I|5b(l . in terms of observable 
quantities, namely, the discontinuities of the magnetiza- 
tion and the magnetic susceptibility, respectively, across 
the first-order transition at T = h = 0. The former is 
given by mi, and differentiating the equation of state 
shows the latter to be A% =1/ 4wm 1 2 . We find 

/i c =le- r / 4 m 1 /Ax»0.23m 1 /A X . (7) 

A very rough estimate using the data from Ref . y predicts 
a value of h c on the order of 0.1 T for ZrZn 2 . 

We now turn to the critical behavior at the various 
critical points in the phase diagram. Equation © yields 
mean-field critical behavior for all points on the lines of 
critical points at T > T tc and T tc > T > 0. In particular, 
the order parameter critical exponents (3 and 8 have their 
mean-field values (3 = 1/2 and 8 = 3, respectively. At the 
tricritical point one finds mean-field tricritical behavior 
[HEU, with (3 = 1/4 and 8 = 5. This behavior gets 
modified if fluctuations are taken into account. For T > 
T tc the upper critical dimension dj, above which mean- 
field critical behavior is exact, is d+ = 4. For d = 3 the 
exact critical behavior is in the classical Heisenberg, XY, 
or Ising universality class, depending on the nature of 
the ferromagnet. (ZrZn 2 is a Heisenberg magnet; UGe2 
has a strong spin anisotropy and is thus Ising-like; MnSi 
is a weak helimagnet, which leads to some complications 
which we ignore here 0.) For the wing-critical lines at 
T < Ttc, d+ = 4 as well. The exact critical behavior is 
always in the Ising universality class, since the external 
magnetic field reduces the effective dimension of the order 
parameter to one. At the tricritical point, d+ = 3, and 
the mean-field theory yields the exact critical behavior 
except for logarithmic corrections to scaling 0] . 

For the quantum critical behavior at T = the mean- 
field theory yields the usual mean-field values for the 
static exponents, e.g., /3 = 1/2 and 5 = 3. For the 
temperature dependence of the order parameter at the 
critical point one finds 6<p(t c ,hc,T) oc -T 2 / 3 , where 
Sip = tp — ific- In the light of Ref. ^3 one might ex- 
pect the exact critical behavior to differ strongly from 
these results. However, a detailed analysis shows that 
Hertz theory holds at this quantum critical point. 
The reason is that the nonzero magnetic field and magne- 
tization suppress the soft-mode effects which invalidate 
Hertz theory, and ultimately destroy the quantum criti- 
cal point, at h = 0. More generally, it was shown in Ref. 



Eol that Hertz theory is valid if the field conjugate to the 
order parameter does not change the soft-mode structure 
of the system. In the present case, an expansion in pow- 
ers of Sip about the quantum critical point shows that 
the quantity 2ip c St — 8h, with St = t — t c and Sh = h — h c , 
plays the role of the conjugate field. Switching on an ex- 
ternal magnetic field from h = gives certain soft modes 
a mass, but changing h from h c 7^ does not lead to fur- 
ther changes in the soft-mode spectrum, and neither does 
changing the value of t. Mean-field theory thus gives the 
exact static quantum critical behavior, in particular 

(3=1/2 , 8 = 3. (8) 

However, the dynamic quantum critical behavior, i.e., the 
temperature dependence at criticality, is modified from 
the mean- field result 0] , since the leading temper- 
ature dependence of the parameter t appears only at 
one-loop order. This fluctuation effect leads to a tem- 
perature scale with a scale dimension [T]a uc = 9/(d+l). 
For d < 5 this dominates the Fermi-liquid temperature 
scale, which has [T]fl = 3/2 and is responsible for the 
temperature dependence of the order parameter within 
mean-field theory. In d = 3 we thus have the exact result 

«V(i c A,T)oc-T 4 / 9 . (9) 

Notice that the static order parameter does not depend 
on the critical temperature scale, which determines the 
dynamical critical exponent z proper, 

z = [T] c = 3. (10) 

For d > 2, the critical scale dominates the fluctuation 
scale for all observables that depend on it, e.g., the spe- 
cific heat . Notice that the above results are the 
exact quantum critical behavior. 

We finally discuss the relation between the theory pre- 
sented above and a competing mean-field theory with 
a very different microscopic underpinning. Sandeman 
et al. ^3 have proposed a Stoner model where the 
equation of state is analytic in the order parameter, but 
band structure effects, in particular a double-peak struc- 
ture in the density of states near the Fermi level, lead 
to signs of the coefficients consistent with a first-order 
transition. These authors have shown that this provides 
an explanation, not just for the first-order nature of the 
paramagnet-to-ferromagnet transition, but also for a sec- 
ond, metamagnetic, transition observed in the ferromag- 
netic phase of UGe2, and they have argued that it also 
leads to triplet superconductivity within the ferromag- 
netic state, in agreement with observations on UGe2 and 
URhGc. Band structure calculations for these two ma- 
terials have confirmed that a double-peak structure near 
the Fermi level exists It is interesting to compare 

various features and predictions of these two theories. 
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(1) The Stoner theory relies on detailed band-structure 
effects to explain the first-order nature of the transi- 
tion. The present theory, on the other hand, is based 
on a universal many-body effect, namely, the existence 
of soft particle- hole excitations, which are always present 
in metals. It therefore predicts the first-order transition 
to be a generic feature of low-T c itinerant ferromagnets, 
independent of the details of the band structure. 

(2) Within the Stoner theory one expects a temperature 
dependence of the coefficient u in Eq. @ from Fermi liq- 
uid theory @, namely, u = u - m (T/T ) 2 . Here T 
is the same microscopic temperature scale as in Eq. (Q, 
and uq ju\ is on the order of unity. One therefore expects 
Ttc to be generically on the order of To, and it is a pri- 
ori not clear what suppresses T tc to the observed values 
around 10K. The many-body theory, on the other hand, 
provides a natural explanation for this effect: The coef- 
ficient v in Eq. reflects a mode- mode coupling effect, 
and therefore v/u <C 1 0- T tc is thus exponentially small 
compared to To. 

(3) Both theories yield magnetic-field dependences of the 
phase diagram that are qualitatively the same, and quan- 
titatively very close to one another. For instance, the 
relation given by Eq. (JJJ) is the same in the Stoner the- 
ory, only the coefficient changes to 3 x 2 7 / 2 /5 5 / 2 w 0.61. 
Notice that no magnetic field dependence of the coeffi- 
cients of either theory is necessary in order to produce the 
characteristic "wing structure" of the phase diagram, the 
term —h<f> in the free energy suffices. In fact, the "wing 
structure" is a direct consequence of the existence of the 
tricritical point and will be present in any theory 
that describes the latter. 
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